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We consider the behavior of an ultracold Fermi gas across a narrow Feshbach resonance, where the 
occupation of the closed channel may not be negligible. While the corrections to the single-channel 
formulae associated with the nonzero chemical potential and with particle conservation have been 
considered in the existing literature, there is a further effect which to the best of our knowledge has 
until now escaped attention, namely the "inter-channel Pauli exclusion principle" associated with 
the fact that a single hyperfine species may be common to the two channels. We focus on this effect 
and show that, as intuitively expected, the resulting corrections are of order Ef In, where Ef is the 
Fermi energy of the gas in the absence of interactions and rj is the Zeeman energy difference between 
the two channels. We also consider the related corrections to the fermionic excitation spectrum, 
and briefly discuss the collective modes of the system. 



I. INTRODUCTION a, > o, s < a 



In a low-temperature dilute system, a short-range in- 
teraction can be characterized with a single parameter, 
a s , a.k.a. the s-wave scattering length. A very desirable 
property of the Feshbach resonance in such a system is 
that the effective interaction is tunable experimentally 
through the Zeeman energy difference between channels 
which is in turn tunable through instruments such as a 
magnetic field [D-ll]- 

a s (B) = a hg (l + ^^) (1) 

where Bo is the magnetic field at which the a s diverges, 
i.e., the resonant point. This unique ability gives physi- 
cists a rare opportunity to study a many-body system un- 
der various interaction strengths, and thus connect differ- 
ent physics originally developed separately. Particularly 
for the fermionic gas, there are a series of theoretical 
works about uniform treatment of BEC and BCS since 
the 1960s [5-8], for which the dilute ultracold fermionic 
alkali gas with Feshbach resonances provides a perfect 
testing ground. Indeed, these theories work quite well 
qualitatively. 

The two-body theory of the Feshbach resonance has a 
characteristic parameter, S c ~ a% (AB) 2 , defined as the 
detuning energy at which the weight of the bound state 
shifts from predominantly in the open channel to pre- 
dominantly in the closed channel (see Fig. [IJ [9J . Naively 
speaking, on the negative detuning side of any resonance 
(i.e. S < 0), the two particles should mostly stay in a 
(virtual) bound state of the closed channel (or "virtual 
state" in some other- type resonances). However, at the 
resonance point of a Feshbach resonance (a s — >■ ±oo), the 
atoms are mostly still in the open channel, and they do 
so down to a negative detuning 5 ~ —S c . Only when the 
negative detuning from resonance is much more than S c , 
do atoms have the majority weight in the closed channel. 
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FIG. 1: Energy levels in a Feshbach resonance 
S is the energy detuning from the resonance point. 
The horizontal line stands for the zero energy s-wave 
scattering state, tp ~ - — — , which exists for any 
detuning. The lower solid line stands for the real 
bound state, which only exists for negative detuning 
(S < 0, a s > 0). The dashed line stands for the 
(uncoupled) closed-channel bound state. An 
interesting point to notice is that the real bound state 
appears earlier than the cross point of the (uncoupled) 
closed-channel bound-state level and zero energy. 
Another important point to notice is the negative 
detuning — S c . When the negative detuning is smaller 
than S c , this real bound state is composed mostly with 
atoms in the open channel and vice versa. 



When considering a many-body system with a Fesh- 
bach resonance, an important question is how this en- 
ergy scale, S c , compares to the typical many-body en- 
ergy scale, namely, the Fermi energy of the free fermionic 
atoms, Ep. In the region not too far away from the reso- 
nance (| 5 1 <C S c ), the closed-channel weight is negligible if 
the Fermi energy is much smaller than <5 C , (i.e., broad res- 
onance). Crossover experiments are usually performed at 
detuning not too far from the resonance, and hence the 
closed channel can be safely ignored at the many-body 
level. Eventually, when the detuning is too far away, 
| S | » S c , the bound state is almost like the uncoupled 
closed-channel bound state with a little dressing from 
the open channel. Nevertheless, such a large detuning 
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in the broad resonance is not very interesting because 
the resonance effect is very small and the s-wave scat- 
tering length, a s , is close to its background value then. 
Given the above consideration, for most purposes, we can 
almost neglect the closed channel and the problem can 
be well-described as a two-species fermion system with a 
tunable interaction when it is not too far away from the 
resonance. The Feshbach resonance indeed serves as a 
simple "magic" knob to change the interaction strength. 
The original theories developed on single-channel models 
apply to this case directly. This is also the situation for 
two popular experimental systems ( 6 Li atoms at 834G, 
40 K atoms at 224G). Many theoretical works have been 
developed using either the single- channel model or the 
two-channel model with the broad resonance assumption 
(e.g. [lol - [l2l |). On the contrary, when the Fermi energy of 
the free fermionic atoms is comparable to or even larger 
than <5 C , the closed channel has to be included at the 
many-body level even for small detuning. Such a situ- 
ation, previously considered in some works [13j, is the 
focus of the current work. 

To complicate the problem further, configurations of 
Feshbach resonances often have one common hyperfine 
species between the two channels. There are three hy- 
perfine species in the two channels instead of four species 
(two for each channel). Two most common systems ( 6 Li 
at 834G, 40 K at 224G) both contain three species of 
fermions although they are broad resonances. The Pauli 
exclusion principle prevents atoms of both channels from 
occupying the same momentum level simultaneously be- 
cause of this common species. This "inter-channel Pauli 
exclusion" has no counterpart in the two-body physics. 
This peculiar effect in many-body crossover problems has 
received little theoretical attention^ up to now. Neverthe- 
less, narrow resonances do exist [3] and it is not incon- 
ceivable to perform many-body experiments using such 
resonances. The central concern of this paper is about 
these situations. 

Roughly speaking, turning from two-body systems to 
many-body systems brings three new effects into the orig- 
inal two-body problem. The first effect is closely associ- 
ated with the Fermi energy: For a many-body fermionic 
system at low temperature, most fermions are inactive; 
only the fermions close to the Fermi surface participate 
in the interaction processes. Therefore, the energy often 
needs to be measured from the Fermi surface instead of 
from zero as in a two-body situation. The second effect 
relates to counting. Unlike in the single- channel prob- 
lem, there are two relevant densities in the two-channel 
problem: the density of atoms in the open channel, n D , 
and the density of atoms in the closed channel, n c . When 
the closed-channel weight is small (broad resonance), it 
is legitimate to treat the total density as the same as 
the open-channel density. However, in the narrow reso- 
nance, where the closed-channel weight is not negligible, 
counting becomes complicated. Extra care is required 
to specify which channel quantities such as "density" re- 
fer to. These two aspects have been extensively studied 



previously [Tj|. 

The last effect is unique to the three-species problem, 
where one common species is shared by both channels. 
The phase spaces of two channels overlap because of the 
common species, which prevents both channels have the 
occupation in the same momentum level simultaneously. 
This effect is controlled by the wave-function overlap of 
the states in the two channels. A rough estimate of 
this overlap can be made: The uncoupled closed-channel 
bound state which is in resonance with the open-channel 
zero energy threshold has relatively small spatial exten- 
sion, a c . Its binding energy Ef, is close to the Zeeman 
energy difference between two channels, 77. On the other 
hand, fermions in the open channel fill the lowest momen- 
tum states up to typically the Fermi energy, Ep. By a 
simple dimensional argument, the ratio Ep/rj must con- 
trol the overlap effect. This effect has not been addressed 
in any theoretical work to out knowledge. How it mod- 
ifies the many-body picture is the central topic of this 
paper. 

The present paper is divided as follows: Section [H] de- 
fines the many-body model and introduces the appro- 
priate notation for the 3-species case. Section Mil lists 
the mean- field results and renormalization. Section [ill Al 
lists the results for the fermionic excitation, and section 
IIVI discusses the bosonic modes. We conclude and dis- 
cuss our approach in Section [V] Further details of the 
calculations may be found in ref. (14j . 



II. MODEL AND NOTATION 

We denote the three hyperfine species as a, b and c, 
where the open channel contains the pair of species (a, b) 
and the closed channel contains the pair of species (a, c). 
We denote the Grassmann variable for species i by ipi 
and write the two channels in the notation: 



(ipip) = (i>ai>b ipai>c 



i>bi> a 



(2) 



The two-body interaction can then be written as a (2 x 2) 
hermitian matrix U 



U 



U Y 
Y* V 



(3) 



We can now write the finite-temperature action as 



dr / d d r 



2m 



(4) 



Here r\i is the Zeeman energy for hyperfine species i. Wc 
choose the zero so that r\ a = T)b = 0, rj c = V- We per- 
form the Hubbard-Stratonovich transformation here. In- 
troduce 2-component auxiliary fields (order parameters) , 
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(Ai, A2), coupled to the fermionic fields as 

A — ► A - U(i/>ip) 
We can introduce a spinor representation 

AO 



(5) 



(6) 



The action can then be rewritten in a more compact form 
with respect to \I/ and *f> 



The Matsubara summation can be performed by the nor- 
mal trick of multiplying the summand by a Fermi func- 
tion and deforming the contour [la ]. For the summation 
at zero temperature, we just need to consider the positive 
roots, -Eik- It is straightforward to find 



E 



(£ik-6c)(£ik + £k + r?) 



A 2 
^2 



closed 



E 



(Ei k + E 2 k) (Ei k + E 3 k 

(-^lk — £k)(£dk + £k 



A? 



(La k + E 2 k) (-Ei k + E 3 k) 



(13) 
(14) 



S(A,A,V^ 



dr / cf r 



(7) 

where the fermionic correlation Q 1 in the momentum- 
frequency representation is 

f iu n - ik Ai A 2 \ 

Ai + & (8) 

A 2 iUn + ^k+ll) 

where £k = h 2 k 2 /2m — [i. The action in Eq. [7]is bilinear 
in the quantities "J, ^ and we can formally integrate them 
out, with the result 

5(A,A) = /dx (Al^A-trlnCr 1 ) (9) 



III. MEAN-FIELD RESULT AND 
RENORMALIZATION 

Eq. ((SJ can be inverted to get G. The final mean- 
field equations are (for simplicity, both Aj's are taken as 
real. }15l|) 



Ai 
A 2 



U Y 
Y* V 



E 



'2k 



(10) 



where hik and /i 2 k are expectations of the abnormal 
Green's function. 



hiw = (V>a,-k^6,+k) = Ai 



fok = (Vfe.-k^c.+k) = A 2 



Eik + £k + r? 



(Ei k + E 2 k) (Ei k + E 3 k) 



(11) 



El k + £k 



(Ei k + E 2 k) (Ei k + E 3 k) 



(12) 



where L^k's are the eigenvalues of the fermionic correla- 
tion Eq. E| (see details in Sec. IIIIA|) . 

There is one number equation for each channel, 



J2 G 22 e ( - ,; "" 5 - ) = A", 



open 



G 33 e(- iu ^ = N close 



A. The Bogoliubov canonical transformation and 
the fermionic excitation 

The fermionic correlation (Eq. [5]) can be diagonalizcd 
with a (unitary) Bogoliubov canonical transformation. 
We break down the transformation into two steps Tk and 
L k . 



^„,k = 4^G- 1 



k^kEk 



(15) 



Here Lk is the diagonal matrix; Tk and Lk are both 
unitary transformations. We take Tk as the canonical 
transformation at the broad resonance, i.e., when we can 
ignore the inter-channel Pauli exclusion. 



u k v k N 
-v k u k 
1, 



(16) 



where u k and v k are defined in a similar fashion as in the 
single-channel BCS problem 



1 



2 = l( 1 _6 L 

Lk^^ + A 2 ) 1 / 2 



(17) 
(18) 



Note that in the narrow resonance here, v 2 . does not carry 
the physical meaning of the occupation number of the 
(open-channel) atoms, and Lk does not represent the 
fermionic excitation spectrum. Tk can be taken as the di- 
rect sum of two parts, the first two rows/columns describe 
the open-channel excitations, while the third row/column 
describes the closed-channel ones. It nevertheless is not 
sufficient to diagonalize the fermionic correlation in the 
narrow resonance. 

Introduce a dimensionless scale £, 



c = 



A 2 

A177 



El 

E b 



1/2 



< 1 



(19) 



where L& is the binding energy of the uncoupled closed- 
channel bound state, cf>. Here both Ai and A 2 are their 
mean-field (saddle point) values. It is not hard to find 
the additional unitary transformation Lk to first order in 
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/ n AjA 2 




A, 



= I+S k (20) 



So we finally arrive at the diagonal matrix B Unt k to first 
order in £ 



B Un ,k = iu n I - ( 





(21) 



The eigenvalues of Q 1 in B Unt k describes the dispersion 
spectrum of the fermionic excitations 



E2W 
Esk 



Ek 



Ek~ 



& + »? 

A 2 ,, 2 
^2% 



2(£k + v) 



Ek + KAiC 



Ek - «kAiC 



£k + V 



(22) 
(23) 
(24) 



i?ik and i?2k correspond to the traditional Bogoliubov 
quasi-particle modes; while E^k describes the fermionic 
excitation mostly in the closed channel. The correction 
due to the inter-channel Pauli exclusion are all of order 

C- 



B. Renormalization of the mean-field equation 

The summations in the gap equations above (Eq. 1 10[) 
diverge when they are converted into integrals at 3D be- 
cause of the artificial assumption of contact interactions. 
We can nevertheless remove the singularity in two steps. 
First, we notice that the closed-channel bound state is 
much smaller than the interparticle distance. Therefore 
the two-body correlation within the closed channel is al- 
most unchanged from its two-body value. Given this 
consideration, we project the closed-channel correlation 
/12k onto the two-body bound-state wave function of the 
uncoupled closed channel, 4>. Second, the bare interac- 
tion in the open channel is replaced with a more physi- 
cally meaningful quantity, the effective s-wave scattering 
length, a s , which already incorporates the singularity. 
We refer to ref. [14[ for the details of the calculation. 

First, we project the closed-channel correlation /12k 
into cbn 



'2k 



ac 



'Ok u k 



(25) 



This projection is only for the high momentum, for the 
low momentum, the available phase space for the closed 
channel is much more limited because the atoms in the 
open channel center in the low momentum. More specif- 
ically, this restriction is represented by the factor 



We come to this conclusion because the closed-channel 
bound state <f> is much smaller in size than the interpar- 
ticle distance. This guarantees that the low momentum 
states are dominated by the open channel. Comparing 
this equation with the two-body Schrodinger equation for 
d>, we find 



E 



kk' 



^k^kk'/llk' 



(-£ b + 77 - 2 M - Ai) 



(26) 



where 



X 1 (r 1 ) = -J2K(Ek- 6<)0k - <^k<Vkk^k' (27) 



kk' 



Put all these together into the gap equation of the open 
channel, we find 



Aip=E( C7 > 



'pk 



EkV Y pp >(t) p >(f>* k ,Ykk> 
(-E b + r]-2fi- Ai) 

Ekk'p' Y pp ,d> p <djl,Ykk>vl,hik 



fin 



(-E b + r ] -2 f i-\ 1 ) 

Comparing to the similar equations in the two-body 
problem, we find the detuning from the resonance here 
differs from that of the two-body physics by 2/i + Ai. 2fi 
describes the many-body shift of the starting point from 
zero to the Fermi surface, while Ai describes the inter- 
channel Pauli exclusion. Furthermore, the last term in 
the r.h.s. is also unique for the many-body problem. This 
term has no singularity for high-momentum because of 
the extra v p , term in the integral. 

Now we handle the above equation with the similar 
strategy as in the single- channel problem: we integrate 
out the high-momentum part and replace the bare in- 
teraction with the more physically observable effective 
open-channel a s . Multiply both side with (1 + TG), 
where T is the scattering matrix for the open channel, 
and G = (u> — Hq)^ 1 is the Green's function for a free 
pair in the open channel. fll\ 



(l + TGQAi = -Tfci-Aa 
A2 P = + TG) pp Y pp/ (f) p/ i 



p' 



(28) 



pp 



Using the zero energy value of the free pair Green's 
function G(lu = 0) = (— 2ek) _1 and introducing the low- 
momentum s-wave scattering length with extra shift 

a s =a bg (l + ^^) (29) 

where /C = A M S, A M is the effective difference of the pair 
magnetic moments of two channels. And 5 is the energy 
detuning between two channels, we have the renormalized 
open-channel gap equation 

Ax 



1 



' [ 2Ek 2e k 



4(& + V)K 



(30) 
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Here the two factors Ai and A2, given respectively by 
Eqs. (|27l) and (|28p. as well as the last term in the 
above equation, describe the inter-channel Pauli exclu- 
sion effect. All of them involve overlap integrals between 
the open-channel wave function and the closed-channel 
wave function. (The factor d£ or w k describes mostly 
the open-channel wave function, while <f> describes the 
closed-channel wave function). The larger the overlap 
of the two, the larger are Ai and A2. This has a very 
intuitive interpretation: more overlap leads to more se- 
vere inter-channel Pauli exclusion, which in turn leads 
to larger (corrections) terms. In our model, the open- 
channel wave function is spread over a large region of 
real space, (even on the BEC side, the real bound-state 
is very loosely bound comparing to the closed-channel 
bound state), while the closed-channel wave function is 
very sensitive to the binding energy, Bfe(~ rf). A closed- 
channel bound state is more spread out in real space and 
has larger overlap with the open-channel wave function, 
if it is closer to the threshold, i.e., the binding energy is 
smaller. Consequently, the terms Ai and A2 are larger 
in such cases. Nevertheless, Ai is much smaller than the 
Fermi energy Ep, or the other shift, the chemical poten- 
tial, 2(i. So it is still legitimate to treat this shift as a 
perturbation. In addition, it is not hard to see that all 
these terms are linear in the density. 

We can also use the expansion on E^ in Eqs. (I22j|24[) 
to rewrite the two number equations Eqs. (|131 1141) to 
first order in £ 



2fi in a s (Eq. 



N, 



closed 



E 



A 2 



Kk + T?)(2& + 7/) 



N, 



k 



2E k V 77 w 4£3\ 



C. Discussion of the mean-field solution 



As discussed before, the correction of the narrow Fes- 
hbach resonance can be taken into account in two steps. 
First, omitting the inter-channel Pauli exclusion, we only 
consider the chemical potential fi in the shift and the ex- 
tra counting due to the closed channel. Then in the sec- 
ond step, we can correct the previous result with quan- 
tities originated from the inter-channel Pauli-exclusion 
unique to the three-species problem. 

In the first step, the gap equation and the (open- 
channel) number equation are simplified to 



1 = - 



4ira s (p) 



E(- 



E 



1 

2-Ek 2ek 
Ek — 6c 
2E V 



(33) 
(34) 



a bg (l 



K 



6-2\i' 6-2fi 



(35) 



The above equations need to be solved self-consistently. 

An interesting point to notice is that the gap Ai sat- 
urates in the BEC side. This is because the effective at- 
traction in the open channel where the pairing happens 
is limited by the real attractive strength in the closed 
channel. It can not become infinitely strong as in an 
ideal single-channel model. Mathematically, this can be 
seen from the gap equation, Eqs. (|10H12[) . The closed- 
channel correlation is limited by the total density, and 
therefore the gap has a maximum. 

Once this step is finished, we can look for the correction 
due to the inter-channel Pauli exclusion. Both Ai and A2 
can be shown to be linear with the density, so can be 
calibrated by experiments of different densities. With 
the proper value of these two parameters, we can find 
the correction numerically. It is nevertheless not hard 
to show that the correction is of order (, so that we are 
warranted to treat this effect MS M perturbation. 



IV. BEYOND MEAN-FIELD: COLLECTIVE 
MODES 



The order parameters (Ai, A2) are defined in terms 
of the collective behaviors of many fermion atoms. Fluc- 
tuations of the order parameters thus signal the collec- 

(31) 

tive excitation modes of the system. Here with a two- 
component order parameter, four independent modes ex- 
ist: two for the magnitude variation of each Aj, the in- 
ternal phase between two A, , and the overall local phase 
(32) ®( x ) °f Ai and A2. The first three change the magni- 
tude of the action and therefore are massive; while the 
last one leaves the action invariant and thus is massless. 
We proceed to study two modes of the phase fluctuation. 

The in-sync phase mode is the counterpart of the 
Anderson-Bogoliubov modes in the single-channel prob- 
lem. The action 5(Aj,Aj) (Eq. [9J, is invariant if the 
phases of A^ and A 2 k rotate simultaneously. We there- 
fore conclude that there exists a massless (Goldstone) 
mode corresponding to the local phase invariance. Intro- 
duce the phase fluctuation 9, 



Ai(x) ->• A 4 e 



.26 (x 



-i20(a 



Here we only consider the shift of the chemical potential 



We have a new action for the phase fluctuation 8 

S[9] = j dxY j e{q)e{-q)[\-K^{Q)wl ~ £- 9 2 ] (36) 

The fluctuation has a linear dispersion relation and there- 
fore is a sound-like mode. It is not hard to show that the 
velocity follows the same structure as that for the single- 
channel crossover problem Q with a correction of order 

C- 
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The out-of-sync phase mode can be studied with the 
similar approach. However, for this mode, we can show 
that the action is 

s[8] = J dx^oim-q^Ho)^-^) £-, 2 ] 

(37) 

where uj m has a finite value loq at the zero momentum, 
which indicates a gapped mode. It can be shown that ujq 
is in the order of the ionization (pair breaking) energy, 
which is around Ai on the BCS side, while around on 
the BEC side. 



V. CONCLUSION AND DISCUSSION 

In this paper, we have studied the narrow Feshbach 
resonance in the three-species case where two channels 
share the same species. 

In general, for the narrow resonance without a shared 
species, the main correction to the single-channel result 
comes from the extra counting of the atoms in the open 
channel, which leads to the extra shift 2/i in a s , and the 
closed channel, which leads to the extra number equa- 
tion. Two number equations exist, one for each channel. 
The open-channel number equation resembles the num- 
ber equation of the single- channel model. 

When there is a common species, however, the Pauli 
exclusion between two channels due to the common 
species in the three-species narrow resonance, calls for 
careful consideration. 

Our treatment follows the the idea of "universality" 
[l~8l HU . For a dilute system with a short-range poten- 
tial, such as the dilute ultracold alkali gas, the short- 
range part of a two-body correlation does not signifi- 
cantly change from two-body to many-body. This par- 
ticular feature justifies using the two-body quantities as 
the boundary condition for the many-body correlations. 
In non-resonant situations, the s-wave scattering length, 
a s , is a good candidate for such a purpose. However, in a 
Feshbach resonance, the open-channel scattering length 
becomes enormous, comparable to or even larger than 
the many-body length scale; furthermore, its value is very 
sensitive to the ratio of weights of atoms in the two chan- 
nels. Therefore it is no longer suitable as the two-body 
characteristic quantity used in the many-body problem. 
A better candidate of two-body quantities as boundary 
condition here is the short-range part of the two-body 



wave function itself. In other word, we expect the two- 
body correlation of the many-body system matches the 
shape of the two-body wave function at short distances. 

When the spatial extension a c of the closed-channel 
bound state is of the order of the inter-particle distance, 
ao, or even larger, the Feshbach resonance in the many- 
body context is a genuine three-species many-body prob- 
lem and no simple solution is available. By contrast, 
when the bound-state's spatial extension is much smaller 
than the interparticle distance, the two-body correlation 
in the closed channel, which is expected to almost entirely 
concentrated in the short-range part, is proportional to 
the two-body bound-state wave function. The ratio of the 
the bound-state size and the interparticle distance, a c /ao, 
serves as the expansion parameter and we can extract the 
effect of the inter-channel Pauli exclusion perturbatively. 
In essence, we can then ignore the many-body effects 
within closed-channel bound states, while only taking 
into consideration the Pauli exclusion between channels 
and within the open channels. A few new parameters 
need to be introduced and can be calibrated from exper- 
iments, such as Ai, A2. Mean field properties can still 
be determined through gap equations and number equa- 
tions similarly to the single- channel case. The excitation 
modes are also close to the original single-channel result 
with correction of the order of £. 

In our approach, we take the broad resonance result 
(or the single-channel crossover) as our zeroth order so- 
lution, upon which the expansion is performed. It is 
however known that the simple BCS-type pairing treat- 
ment is not adequate to quantitatively describe the whole 
BEC-BCS crossover region. Therefore the zeroth order 
solution used here can be improved through further theo- 
retical development. Nevertheless, we expect the pertur- 
bation approach used here to build the narrow resonance 
from the single-channel crossover result to remain valid. 
Once the zeroth order solution (for a broad resonance or 
a single channel BEC-BCS crossover model) is appropri- 
ately improved, the correction of the narrow resonance 
in such a parameter regime, can still be obtained by a 
procedure similar to the one in this paper. 

ACKNOWLEDGMENTS 

We thank Professor Monique Combescot, Dr. 
Shizhong Zhang and Dr. Wei-Cheng Lee for many in- 
spiring discussions. Part of this research is supported by 
the National Science Foundation under grant No. DMR 
09-06921. 



[1] U. Fano, |Phys. Rev. 124, 1866 (1961)| 
[2] H. Feshbach, Annals of Physics 19, 287 (1962). 
[3] C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, 
|Rev. Mod. Phys. 82, 1225 (2010) | 



[4] C. J. Pethick and H. Smith, Bose-Einstein Condensation 

in Dilute Gases (Cambridge University Press, 2001). 
[5] D. M. Eagles, [Phys. Rev. 186, 456 (1969)| 
[6] A. J. Leggett, in Proceedings of the XVIth Karpacz 



7 



Winter School of Theoretical Physics, Karpacz, Poland 

(Springer- Verlag, 1980) pp. 13-27. 
[7] P. Nozieres and S. Schmitt-Rink, 

|J. of Low T. Phys. 59, 195 (1985)] 

10.1007/BF00683774. 
[8] M. Randeria, in Bose-Einstein Condensation, edited by 

A. Griffin, D. Snoke, and S. Stringari (Cambridge Uni- 
versity Press, 1995) pp. 355-392. 
[9] A. J. Leggett, Quantum Liquids (Oxford University 

Press, 2006) appendix 4A. 
[10] M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chio falo, 

and R. Walser, |Phys. Rev. Lett. 87, 120406 (2001) | 
[11] T.-L. Ho, |Phys. Rev. Lett. 92, 090402 (2004) 
[12] J. N. Fuchs, A. Recati, and W. Zwerger, 

|Phys. Rev. Lett. 93, 090408 (200!) 
[13] V. Gurarie and L. Radzihovsky, 

|Annals of Physics 322, 2~ ( 2007) , january Special 

Issue 2007. 

[14] G. Zhu, BEC-BCS Crossover with Feshbach Resonance 
for Three- Hyperfine- Species Model, Ph.D. thesis, Univer- 
sity of Illinois at Urbana-Champaign (2012). 

[15] When Y is not real, Ai and A2 cannot be both real even 



at the mean field level. Nevertheless, we can require one 
real, then the other will have a phase just to compensate 
the phase in Y . The final conclusion can be verified to 
remain valid. 
[16] See sec. 4.2.1 in 0, sec. 25 in HH- 

[17] Here we use the relation between the scattering T-matrix, 
the free pair Green's function G — (uj — Ho)" 1 and the 
bare interaction V = — U e g. 

T = V + TGV = -(1 + TG)U cff 

[18] S. Tan, [Annals of Physics 323, 2952 (2008)] 

[19] S. Zhang and A. J. Leggett, 

|Phys. Rev. A 79, 023601 (2009) 
[20] A. Altland and B. Simons, 

Condensed Matter Field Theory (Cambridge University 

Press, 2010). 

[21] A. L. Fetter and J. D. Walecka, Quantum Theory of 
Many-Particle Systems, dover ed. (Dover Publications, 
Inc. (Original Publisher: McGraw-Hill), 1971). 



